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An ultracold gas of heteronuclear alkali dimer molecules with hyperfine structure loaded into 
a one-dimensional optical lattice is investigated. The Hyperfine Molecular Hubbard Hamiltonian 
(HMHH), an effective low-energy lattice Hamiltonian, is derived from first principles. The large 
permanent electric dipole moment of these molecules gives rise to long range dipole-dipole forces 
in a DC electric field and allows for transitions between rotational states in an AC microwave 
field. Additionally, a strong magnetic field can be used to control the hyperfine degrees of freedom 
independently of the rotational degrees of freedom. By tuning the angle between the DC electric 
and magnetic fields and the strength of the AC field it is possible to control the number of internal 
states involved in the dynamics as well as the degree of correlation between the spatial and internal 
degrees of freedom. The HMHH's unique features have direct experimental consequences such as 
quantum dephasing, tunable complexity, and the dependence of the phase diagram on the molecular 
state. 



I. INTRODUCTION 

Ultracold molecular gases are of interest in many sub- 
fields of science ranging from precision science to quan- 
tum simulation of many-body Hamiltonians [ij. Recent 
success using the STIRAP (STImulated Raman Adia- 
batic P assa ge) method has allowed experimentalists to 
produce a gas of KRb molecules close to Fermi degener- 
acy, in the ground rovibrational state, and in a specific 
hyperfine level [3, ^ . Rovibonic ground state molecules 
have also been formed for polar LiCs[^] as well as nonpo- 



lar Cs2[5| and Rb2 
rently underway [3, 



with studies on other species cur- 
To reach the quantum degenerate 
regime one must have all molecules in the same quantum 
state, a task which is complicated by the rich hyperfine 
structure of alkali dimer molecules. Thus, a number of 
recent works (ol-ITTI have investigated the single-molecule 
microwave spectra to find a route by which all molecules 
are transferred to the lowest hyperfine state, yielding a 
gas of absolute ground state molecules. 

From the condensed matter perspective, ultracold 
gases are enticing in their capacity to act as quantum 
simulators [l2|, Ell . Such specialized quantum comput- 
ers allow for the study of complex many-body Hamilto- 
nians in a setting where many parameters are amenable 
to experimental control. From this point of view, it is 
natural to ask how the various degrees of freedom in the 
quantum simulator may be controlled and used as re- 
sources. Theoretical proposals for many-body physics 
using ultracold molecules have so far focused only on the 
rotational degrees of freedom in ^ E molecules with exter- 
nal fields (UIHI or on the hyperfine degree of freedom 
in molecules without external fields [16|- In this work 
we study molecules in strong fields including the ef- 
fects of hyperfine structure and discuss how the hyperfine 
degrees of freedom may be controllably accessed and ma- 
nipulated as a resource for generating complex quantum 
dynamics. 

For molecules it has been shown that the interaction 
of the rotational degrees of freedom with external elec- 
tric fields allows for the tuning of the strength and range 



of the two-molecule interaction potential [1^. Many of 
these results also hold for molecules with hyperfine struc- 
ture, as the rotational and nuclear spin degrees are only 
weakly coupled in strong fields. In particular, the ap- 
plication of a DC field and an optical trapping potential 
gives rise to a purely repulsive dipole-dipole interaction 
between molecules in reduced geometries. Also, it has 
been shown that the combination of a strong uniform 
magnetic field and a suitably chosen microwave field al- 
lows for transitions between particular hyperfine single- 
molecule states, and that this may be used to transfer a 
collection of molecules that have been cooled to the rovi- 
brational ground state but an excited hyperfine state to 
their hyperfine ground state (ol-llH . This idea also works 
in reverse: one can select the states which are involved in 
many-body dynamics with the ground state by judicious 
choice of the field strengths and geometries. The HMHH 
refiects this fact; not only the parameters of the Hamil- 
tonian but also the dimensionality and character of the 
basis are suited to experimental control. 



This article is organized as follows. In Sec. [Ill we intro- 
duce the HMHH, define its parameters, and discuss its 
novel experimental consequences. This section contains 
the main results of the paper. In Sec. Ill CI we derive the 
HMHH from first principles and state the key assump- 
tions underlying its derivation. Finally, in Sec. |IVl we 
conclude. Some details concerning the single molecule 
physics are provided in the appendices in the interest of 
completeness. 
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II. STATEMENT OF THE HAMILTONIAN AND 
EXPERIMENTAL CONSEQUENCES 

The Hyperfine Molecular Hubbard Hamiltonian is 



(T i a ' ' 



(1) 



where dia- destroys a bosonic or fermionic molecule in 
state \a) on the i^^ lattice site, and the bracket notation 
(...) denotes that the sum is taken over nearest neigh- 
bors. The single- molecule basis {|cr)} takes into account 
the hyperfine interactions (Appendix [A]) and static fields 
(Appendix [B]) and the quantum number <j is a composite 
index referring to both rotational and nuclear spin de- 
grees of freedom. The properties and dimensionality of 
this basis can be modified by the geometry and strength 
of the external fields, as will be discussed in more detail 
below. 

The first term in the HMHH represents the energy off- 
set of a molecule in state \a) from a reference ground 
state. The second term describes the tunneling of 
molecules between lattice sites and depends on the ro- 
tational state. The third term describes resonant dipole- 
dipole interactions between molecules on neighboring 
sites. The final term corresponds to transitions driven 
between states |cr) and |<j') by an AC microwave field. 
Here the transition dipole moment between two states 
\a') is dcra' = {o'\di\(j')^ where = d • ei is the 
projection of the dipole operator along the space-fixed 
spherical basis direction ei = — (e^; + iSy) /\/2. 

For ^^K^^Rb, which is the most experimentally rele- 
vant species, the energy scales of the various terms are 
summarized in Table HI The detunings A^r are deter- 
mined chiefly by the linear Zeeman effect, and so are 
tunable by the DC magnetic field, and will be similar 
for other molecular species. The tunneling energy scale 
to- is set by the recoil energy, and so will be similar for 
other alkali dimers. The dipole-dipole energy scale U^a' 
is fixed by the permanent dipole moment, and so will 
change with the molecular species. For example, LiCs 
has a dipole moment roughly 10 times larger than that 
of KRb, and so Uaa' will be of order 25kHz. The scale 
of the AC term is determined by the power of the mi- 
crowave field Eac^ which is readily tunable. The range 
of energies we have quoted represents the most interest- 
ing regime where the basic assumptions of our derivation 
hold. 

In the following sections we will justify the HMHH and 
list the essential assumptions underlying its derivation, 
but we first pause to note some of its novel properties. 



Term 


Energy scale 


A. 


1 — lOOkHz (depends on static field strengths) 


ta 


- 1 kHz 




-250 Hz 


daa'E AC 


- 1 - 50kHz 



TABLE I: Table of energy scales of the Hyperfine Molecular 
Hubbard Hamiltonian. From top to bottom: energy /S.^ of 
internal state |cr), relative to the ground state; tunneling ta] 
dipole-dipole interaction U^j^'] transition dipole moment d^^/ 
due to the AC electric drive ^ac- 



A. Quantum Dephasing 

The first property, which we call quantum dephas- 
ing^ was investigated previously for a molecular Hub- 
bard Hamiltonian involving only rotational degrees of 
freedom [l5|. The effect, which is purely many-body in 
nature, may be summarized in this context as the de- 
struction of coherent Rabi flopping due to the popula- 
tion of many spatial degrees of freedom in a many-body 
system driven at a single-molecule resonance. This effect 
is also of interest in the more general context of oscilla- 
tions in a many-body system that are damped by some 
intrinsic mechanism following a quench p!?!, Il8|. 

Dephasing is strongest when the Rabi frequency is 
on the order of the tunneling energies and the differ- 
ence in tunneling energies for the two internal modes 
is also comparable to these two scales. For a system 
with two single-particle levels and 1 and tunneling en- 
ergies to and ti, respectively, this gives the condition 
to ti \tQ — ti\^ wMch cau be achieved with the 
HMHH for reasonable parameter values. The Rabi oscil- 
lations between the two internal states connected by the 
single molecule resonance damp out exponentially in time 
with an emergent time scale r which can be measured ex- 
perimentally, see Fig. [H Dephasing can be observed in 
the structure factors 



s: 



■1)* ' {Uiaflja') ■ 



(2) 



where L is the number of lattice sites; S^"^ can be mea- 
sured in scattering experiments [l9[ . 



B. Internal State Dependence of Phase Diagram 

The dependence of the tunneling energy to- on the in- 
ternal state <j makes the borders of the phase diagram 
shift strongly (e.g. by a factor of 2). This dependence is 
shown exphcitly in Fig. [2l Thus, by preparing a collec- 
tion of molecules in multiple internal states one can study 
interactions of many-body systems in different quantum 
phases and possibly far from equilibrium. Possibilities for 
quantum statics include studying the properties of phase 
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FIG. 1: (Color online) Quantum dephasing in the HMHH. 
The plot shows the behavior of the total number in state 0: 
(no) = (X)^ n^o) when the system evolves under the Hamil- 
tonian ([T]). Quantum dephasing produces an emergent ex- 
ponential envelope on the Rabi oscillation pattern between 
states and 1. Only the number of state is shown for clar- 
ity. The dashed red curve is an exponential envelope fit to 
A/'exp (—t/r) with r = 1441.17ms. The nonexponential be- 
havior near t — 200 is due to the finite size of the lattice. 
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FIG. 2: (Color online) Tunneling matrix elements in a DC 
electric field. Tunneling energies (in kHz) of the = 
(solid blue line) and = 1 (dashed green line) rotational 
states and their difference divided by their arithmetic mean, 
2 (ti — to) / (ti + to), (dash-dotted red line) for KRb in a field 
of lOkV/cm as a function of the effective isotropic lattice 
height Tj = a|Eopt|^ (in recoil energy units). The values of 
the polarizability tensor are taken from Ref. [2Q| ]. 



equilibria as a function of population imbalance and ef- 
fective mass (as determined by the tunneling energy) (21} . 
Also, as the difference in tunneling energy between dif- 
ferent modes depends only on the elements of the molec- 
ular polarizability tensor, measuring the borders of the 
static phase diagram for different internal states also pro- 
vides a novel means to measure this tensor. Possibilities 
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FIG. 3: (Color onhne) Geometry of the HMHH. Counter- 
propagating laser beams along the y and z directions create 
an array of ID tubes, and an additional pair of laser beams 
along X creates a lattice potential. A strong DC field ori- 
ents the dipoles along the direction perpendicular to motion, 
and a magnetic field orients the nuclear spins. An AC field 
of circular space-fixed polarization drives transitions between 
internal levels. 



for quantum dynamics include the study of quench phe- 
nomena for interacting many-body systems in different 
quantum phases. 

C. Tunable Complexity 

A final noteworthy property which was not present in 
the molecular Hubbard Hamiltonians previously studied 
is the possibility of tunable complexity. By complexity we 
mean that the system is comprised of many interacting 
degrees of freedom and displays emergent behavior such 
as the dephasing discussed above. Tunability refers to 
the fact that we may alter the number of internal degrees 
of freedom that are accessed dynamically as well as the 
timescale of their relative interactions. The key point for 
tunability is that the electric and magnetic fields affect 
different degrees of freedom: the electric dipole moment 
and nuclear spins, respectively. We illustrate this con- 
cept, and the corresponding geometries and polarizations 
needed for experiments, in Fig. [3l 

In slightly more detail, tunability is achieved as fol- 
lows. In the presence of an electric field aligned along the 
z direction, dipole moments are induced between states 
having the same nuclear spin projection along the field. 
The introduction of a strong magnetic field defines an 
effective axis of quantization for the nuclear spins while 
leaving the rotational structure unchanged because of the 
strong nuclear Zeeman effect, the weak rotational Zee- 
man effect, and the presence of only weak (quadrupole) 
coupling between the rotational and nuclear spin degrees 
of freedom. In the presence of a strong magnetic field 
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that is not collinear with the electric field it is therefore 
possible to induce dipole moments between states with 
different hyperfine quantum numbers. 

Thus, by changing the relative angle between the elec- 
tric and magnetic fields one can control the number of 
states accessible from a particular state. The power of 
the applied AC field determines the interaction scale and 
the Rabi frequency of these dipole couplings, and the 
strength of the magnetic field determines the energetic 
splittings between states, in turn determining the rela- 
tive rates of internal state population. The HMHH may 
therefore be used as a quantum simulator of a quantum 
complex system where the number and timescale of the 
internal components may be dynamically altered. Precise 
measures of complexity and simulations displaying char- 
acteristic behavior in various regimes will be discussed in 
future work 1221. 



where 



III. DERIVATION OF THE HYPERFINE 
MOLECULAR HUBBARD HAMILTONIAN 



We consider the experimental setup shown schemati- 
cally in Fig. [3l Counter-propagating laser beams along 
the y and z directions create a series of ID optical lat- 
tice "tubes." The intensity of these beams is such that 
the tubes are isolated from one another, and the lattice 
spacing is chosen (e.g. by crossed beams) such that the 
dipole-dipole interaction along y and z is negligible on 
experimental timescales. An additional pair of beams 
creates a lattice potential along the x-direction. The ex- 
perimental techniques required to create this setup have 
been well established for ultracold atoms p3l-[25|. In ad- 
dition to the lattice potential there is a uniform DC elec- 
tric field along the z direction, a uniform magnetic field 
which lies in the xz plane, and an AC microwave field 
propagating in z which is assumed to have circular po- 
larization g = 1 in the space-fixed spherical basis. 

In the lattice is an ultracold quantum degenerate gas 
of heteronuclear molecules characterized by perma- 
nent electric dipole moment rotational constant Bn^ 
rotational angular momentum N and nuclear spins 
Ii and I2. Both nuclear spins are taken to be greater 
than one-half, so that both nuclei have nonzero electric 
quadrupole moments. In second quantization the full 
low-energy Hamiltonian for this setup is 



H 



H\n + + ^AC + ^kin + ^opt 



^/,(r) 



^in =Hrot + ^scal + ^tens + ^r-s + ^quad (4) 

2 

=Bn^^ + cA ' I2 + C3I1 • T . I2 + ^ qN . i, 

2 ^ 

i=l 

2 

= - ^rMArN • B - ^ gifiN (1 - en) h ■ B (5) 



- Edc • d , 
Hac = - Eac • d , 



(6) 
(7) 

-f^opt = — E*pj • a (wopt) • Eopt , (8) 
di • d2 - 3 (di ■ euj (br ■ 



2m 



(9) 



-- dvdv' ij\v) V'^(r') Foodr - r'|) V-(r') ^(r), (3) 



The first line of Eq. (|3]) is comprised of single-molecule 
terms. In order, these are J^in, the Hamiltonian govern- 
ing the internal rotational and nuclear spin degrees of 
freedom; the interaction of the molecule with exter- 
nally applied DC electric and magnetic fields; ^^ac, the 
interaction of the molecule with an AC microwave field; 
f^kin, the kinetic energy of the molecule; and i^opt, the 
interaction of the molecule with the optical lattice po- 
tential. The second line of Eq. ([3]) is the two-molecule 
resonant dipole-dipole force. The main assumptions un- 
derlying this Hamiltonian and our subsequent analysis 
are the following. 

First, the individual molecules are assumed to be in 
their electronic and vibrational ground states, and it is 
assumed that none of these degrees of freedom can be 
excited at the large intermolecular separations and low 
temperatures/relative energies that we consider. 

Second, the characteristic trapping potential length is 
chosen large enough compared to the internuclear axis to 
assume spherical symmetry, i.e. a locally constant poten- 
tial. 

Third, we consider only the lowest two rotational lev- 
els. All AC fields will be sufficiently weak to allow this 
assumption. We also work in the rotating wave approxi- 
mation, which requires that the detuning be small com- 
pared to the driving frequency. 

Fourth, we consider all molecules to be in the lowest 
Bloch band. The AC Rabi frequencies are chosen to be 
small (~l-50kHz) in comparison with the lattice band- 
width {^1{)Er ^lOOkHz) to ensure this assumption. 

Fifth, we work in the "hard-core" limit where at 
most one molecule is allowed per site. This is enforced 
by strongly repulsive dipole-dipole interactions on-site, 
caused by our 2:-alignment of the electric field, as sketched 
in Fig [3l We consider the lattice spacing large enough to 
include only nearest-neighbor dipole-dipole interactions. 
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We neglect the effects of chemical reactions or hyperfine 
changing collisions which occur at very short range. 

Sixth, we neglect dipole-dipole interactions between 
molecules in different ID "tubes." For a consistent level 
of approximation this requires the tubes to be separated 
by twice the lattice spacing. This can be achieved in prin- 
ciple using crossed beams to create larger lattice spacings. 

Seventh, we consider only pairwise interactions of the 
molecules, neglecting three and higher-body interactions. 
This is valid for KRb because the permanent dipole mo- 
ment d = 0.566D is rather small. For molecules such as 
LiCs with larger permanent dipole moments, the three- 
body interaction can play a significant role |26| . 

To derive a Hamiltonian of Hubbard type from Eq. (|3]) 
we follow the standard prescription [l3| of expanding the 
field operators of our second-quantized Hamiltonian in 
a Wannier basis of single-molecule states centered at a 
particular discrete position r^: 



(10) 



where z is a site index and a an index denoting the in- 
ternal state of the molecule. The Wannier basis we use 
is the basis which diagonalizes the internal plus static 
field Hamiltonians Hin + and in which all states with 
= 1 rotate with frequency a;, where uj is the frequency 
of the applied AC electric field. With the field operator 
written in this manner, we find the Hubbard parameters 

ta = - Jdrwl{r -Ti) Hkin^Hopt (r - r^+i) , 

(11) 

= Jdrw^^ (r-Vi) 

and 



i^in + Hp 



(r - r,) , (12) 



-daa'EAC= jdYwlir -Yi)HAcW^' {y -Yi) , (13) 

U^a' = jdYdY' < (r - r,) (r' - r,+i) (14) 
X i^DD (r - r') Wcr (r - r^) (r' - r^+i) . 

The detunings Act are determined by the single- 
molecule spectra, which are well-known [ol, [13 • In the 
interest of the present article's completeness, we have 
included appendices reviewing the basic results and ex- 
plaining them in the context of the present problem. In 
the following sections we discuss the remaining Hubbard 
parameters. 



A. Tunneling Energies 

A key component of the realization of many-body 
Hamiltonians using ultracold molecules is the presence 
of a far off-resonant optical lattice which confines the 



molecules in a reduced geometry. The Hamiltonian of 
this interaction is 

^opt = -Ejpt (r, cjopt) • (cjopt) • Eopt (r, cJopt) , (15) 

where Eopt (r, cJopt) is the optical lattice field and a (cc;opt) 
is the polarizability tensor operator of the molecule, eval- 
uated at the optical lattice frequency c^opt- In our no- 
tation, the circumflex accent (the 'hat') denotes an op- 
erator, the tilde denotes a rank 2 tensor, and boldface 
denotes a rank 1 tensor, or vector. This optical poten- 
tial couples to the electronic degrees of freedom and is 
detuned from resonance by an amount several orders of 
magnitude larger than any hyperfine splittings. Thus 
dependence on the hyperfine quantum numbers in negli- 
gible. For tight optical traps, the optical trap potential 
at each well is close to that of a harmonic trap plus a 
small state-dependent tensor shift of the trap frequency 
affecting levels with > due to the polarizability 
anisotropy [15]. 

When the optical potential is combined with the ki- 
netic portion of the Hamiltonian and evaluated in the 
Wannier basis one obtains the tunneling energies. As the 
tunneling energies are independent of the hyperfine quan- 
tum numbers, we can use results obtained in the case of 
only rotational degrees of freedom, derived in our earlier 
work [l5|. Then the tunneling energies in the eigenbasis 
of ^rot, \NMn), are given by 



^NMn 



exp 



-C 



(16) 



where A = 1.397, B = 1.051, and C = 2.121 are fit 
parameters [28], En the recoil energy, and 



IE, 



opt I 



2Aa N{N^1) 



3 (2Ar - 1) (2A^ + 3) 



(17) 



is the effective lattice height for the \NMn) level. Here 
a is the average polarizability and Aa the polarizability 
anisotropy. 

In the presence of a DC field the rotational levels be- 
come mixed, leading to new effective tunneling energies 
which we denote as InMn^ with N and Mn the cor- 
responding zero field values. This hybridization of ro- 
tational levels in principle also allows tunneling events 
which change the rotational state of the molecule, but we 
can ignore such events because the rotational level sep- 
aration is much larger than the tunneling energies. The 
effective tunneling for the TV = and TV = 1, Mn = ±1 
levels is shown in Fig. [2l The scale is set by the recoil 
energy, which is 27r x 1.44 kHz for KRb in a 1054 nm 
optical lattice. 



B. Two-Molecule Interactions 

Heteronuclear molecules posses permanent dipole 
moments, and thus interact via a dipole-dipole inter ac- 
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tion 



where 



ii^DD (R) 



di • d2 - 3 (^di • gr^ (^br • d2) 



(18) 



where R = r2 — ri, e^^ is a unit vector in the direc- 
tion of R, and d^ is the vector dipole operator of the i^^ 
molecule. In the absence of external fields, this interac- 
tion is off-resonant, leading to a van der Waals interaction 
^DD (R) ~ R~^^ but in the presence of electric fields res- 
onant dipoles are induced and the interaction displays a 
resonant behavior in addition to the R~^ behavior. 

The anisotropic nature of the dipole-dipole force has 
been experimentally shown to dominate the rethermal- 
ization behavior of a molecular gas via inelastic colli- 
sions [29|. This is because a "head-to-tail" arrangement 
of molecules leads to an attractive potential, whereas 
"side-to-side" interactions are repulsive. To ensure the 
stability of an ultracold molecular ensemble and to pre- 
vent losses from inelastic collisions it is crucial therefore 
not only to orient the dipoles using a DC field, but also 
to confine the molecules in a reduced geometry. A thor- 
ough discussion of the nature of the two-molecule spectra 
for molecules without hyperfine structure and its im- 
plications for stability in two dimensions is presented in 
Ref. [l^ . Diagonalization of the full two- molecule Hamil- 
tonian is impractical when hyperfine structure is included 
due to the very large matrices that result. Instead, we 
argue based on comparisons of length and energy scales 
that the hyperfine structure is negligible during the col- 
lisional processes which occur in our proposed setup. 

Our reduced geometry is imposed by the optical lattice 
described earlier. Namely, we consider the case where the 
molecules are confined to move only along the x direction 
and a DC field polarized along the z direction orients 
the dipoles such that all collisions are side-to-side and 
repulsive. The dipole-dipole interaction in this geometry 
reduces to 



^ di -\- di ^ d- 



3 d- 



■ di ^ di 



(19) 



where dq = d- eg is the component of the dipole operator 
along the q direction in the space-fixed spherical basis. 
For z-polarized electric field, the only diagonal compo- 
nents are those involving do. The components of the 
interaction involving d±i couple states with AMn = ±1 
that are separated in energy by an amount of order the 
rotational constant for the DC fields we consider (see 
Fig. [5j). Contributions from these components are sup- 

1/3 

pressed at distances greater than tb = {d'^/B) , of 
order a few nanometers. Thus, at the nearest-neighbor 
distance in a 1054nm optical lattice we consider only the 
diagonal elements of the dipole-dipole interaction. This 
restriction gives rise to the two-body term 

^ niaTlja' , (20) 



d„d„> 



(A/2) 



3 • 



(21) 



In Eq. (|2T|) d^ is the resonant dipole moment of state 
\(j) and A is the wavelength of the optical lattice. We as- 
sume that the long-range repulsive diagonal d^ portion of 
the dipole-dipole interaction is strong enough to prevent 
both the occupation of any one lattice site by more than 
one molecule and access to the region where hyperfine- 
changing collisions involving the d±i dipole moments oc- 
cur. 



C. Interactions with static external fields 

The spectral properties of molecules in co- linear DC 
electric and magnetic fields have been elucidated else- 
where in the literature [ol, [lol, [13 , and the basic results 
of the analysis are given in Appendix [B] for the reader's 
convenience. In this section, we focus on the properties 
of such molecules in non-co-linear fields, in particular on 
the dipole moments. 

The behavior of the molecular dipole moments are con- 
trolled by an external DC electric field which mixes ro- 
tational levels of opposite parity and thus orients the 
molecule. However, a DC field does not couple to the nu- 
clear spins. So for a 2:-polarized field the selection rules 
AMi 0, AM2 are enforced, where Mi and M2 are 
the nuclear spin projections along the field direction. In 
contrast, a magnetic field couples strongly to the nuclear 
spins but only weakly to the rotational angular momen- 
tum due to the relative sizes of the ^-factors The 
magnetic field Hamiltonian thus has eigenstates which 
are energetically distinct nuclear spin states with a quan- 
tization axis given by the field direction. It is in this 
sense that we say the magnetic field defines an effective 
axis of quantization for the nuclear spins. Thus, in the 
absence of internal couplings of the rotational and hy- 
perfine degrees of freedom they may be manipulated in- 
dependently: the rotational angular momentum with an 
electric field and the nuclear spin angular momenta with 
a magnetic field. 

The presence of nuclear quadrupole couplings in alkali 
dimer molecules couples states with the same total an- 
gular momentum projection Mp but different rotational 
and nuclear spin projections. For example, in KRb, the 
interaction couples |<j) = = 1, Mn = 0, Mk, ± 1) 
to \a') = \N = 1,Mn = ±l,MK,MRb) with the lat- 
ter accounting for ~ 10% of the state in the absence of 
fields. (The interaction also couples the Mk ± 1 states, 
but the coupling constant {eqQ)j^ is significantly smaller 
than (eqQ)-^^ and so the mixing is negligible in compar- 
ison) [ll[. Clearly, since the TV = state has only one 
projection Mn = 0, the nuclear quadrupole interaction 
leaves this level unaffected. In a DC electric field where 
the rotational levels become mixed, the states correlat- 
ing with the = levels and = 1 levels both display 
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FIG. 4: (Color online) Distribution of dipolar character. The 
colorbar shows the logarithm of the transition dipole moment 
with the ground state, log(g.s.|(ii as a function of the an- 
gle between the magnetic and electric fields Ob and the state 
index(ordered by energy). Changing the angle between the 
electric and magnetic fields breaks the nuclear spin projec- 
tion selection rule and allows for transition dipole moments 
between many states. Only dipole moments greater than 10~^ 
are displayed. 



quadrupole effects, but these effects are still not identi- 
cal. In strong fields the Zeeman effect dominates over the 
quadrupole coupHng, allowing control over the nuclear 
spins that displays a weak dependence on the rotational 
level. 

Thus, a strong magnetic field defines an effective axis 
of quantization for the nuclear spins, resulting in nuclear 
spin states which are superpositions of states in the basis 
with the axis of quantization along the electric field axis. 
This implies that by changing the angle of the magnetic 
field with respect to the electric field, it is possible to 
change the number of states which are coupled by transi- 
tion dipole moments. This is illustrated in Fig. IH which 
shows the logarithm of the transition dipole moment with 
the ground state as a function of the angle between the 
DC and magnetic fields Ob and a state index (ordered 
by energy). The lowest state index denotes the lowest 
energy state in the = 1 manifold. When the fields are 
co-hnear, one state dominates the dipole spectrum. As 
the angle changes the dipolar character becomes spread 
over many states. These transition dipole moments al- 
low the states to couple in an AC microwave field and 
generate complex dynamics. 



D. Interaction with an AC microwave field 

The introduction of an AC microwave field contributes 
to the Hamiltonian in a similar way to a DC field. In ad- 
dition, the inherent time dependence allows for circular 
and linear polarization as well as the possibility of driv- 



ing transitions between internal states. In the absence of 
hyperfine structure, an AC field of spherical polarization 
q couples the \N = 0, Mn = 0) and \N = 1, Mn = q) lev- 
els, leading to an effective two- level system in the Floquet 
picture [ij. In the presence of hyperfine structure, states 
with different total angular momentum projections Mf 
in the > 1 manifolds become mixed due to the elec- 
tric quadrupole interaction. Thus no rigorous selection 
rules can be estabhshed. This complicates the issue of 
addressing single hyperfine states using microwave fields, 
but it also allows the hyperfine state to be changed using 
microwave fields. Addressing a single hyperfine state can 
be achieved by the application of a strong magnetic field 
such as those used in the STIRAP procedure, which de- 
fines the projections sufficiently to suppress transitions 
to non-target hyperfine states ^9i]. In the presence of an 
electric field, this last comment holds only in the case 
where the two fields are co-linear. When the fields are 
not co-linear many states can be accessed from any one 
state via a microwave transition due to the behavior of 
the transition dipole moments in crossed fields, as was 
described in Sec. IIII CI 

We choose the polarization of the AC field to be purely 
circular, ^ac = 1. A component along q = would lead 
to rapid oscillation of the eigenenergies because the do 
moments induced by the electric field couple to the AC 
field, and this complicates the analysis. Furthermore, we 
consider Rabi frequencies which are much less than the 
bandwidth of the optical lattice so that our approxima- 
tion of being in the lowest Bloch band remains valid and 
we are also justified in using a rotating wave approxima- 
tion. The above considerations together with the single- 
molecule AC Hamiltonian 

J^AC = -d • Eac = -4^Ace"'"' + h.c. (22) 
lead directly to the second quantized Hamiltonian 



(23) 



In Eq. ([23|) the label a refers to the eigenstate |<j) of the 
internal plus static field Hamiltonian + H^^ dcra' = 
(<j|(ii |<j'), and E^j is the energy of state |<j). 

Assembling all the many-body terms expressed in this 
basis, we obtain the time-dependent Hamiltonian 

H = YE^Yhi^ -^t^Y^ ^d\^dj^ + h.c. 

cr i a {i^j) 

- i ^ d,,> Eac e'"* + h.c] . (24) 

aa' i 

If we change to a basis where all single-molecule states 
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with = 1 rotate with frequency uj we have, finahy: 
H = ^/^a^ni^ - ^UY^ kl^ja + h.c. 



(25) 



where Act = Ea- for states with TV = and — oo for 
states with = 1. 



IV. CONCLUSIONS 

We have presented and derived the Hyperfine Molec- 
ular Hubbard Hamiltonian (HMHH). The HMHH is a 
lattice Hamiltonian describing the effective low-energy 
physics of an ultracold gas of heteronuclear alkali dimer 
molecules with hyperfine structure loaded into a ID op- 
tical lattice and interacting with external DC electric, 
AC microwave, and static magnetic fields. By tuning 
the angle between the electric and magnetic fields and 
the strength of the magnetic and AC fields it is pos- 
sible to change the number and timescale of internal 
states contributing to many-body dynamics. The Hamil- 
tonian also displays emergent quantum dephasing, and 
has a phase diagram which depends strongly on the ini- 
tial state. These features make the HMHH an ideal can- 
didate for a model quantum complex system. 

Future work will involve time-evolving block decima- 
tion simulations of the HMHH similar to past studies of 
molecular Hubbard Hamiltonians [l5|. In particular, we 
will discuss measures of complexity and how they relate 
to experimentally measurable quantities. Future work 
on the Hamiltonian itself will include realistic models 
of molecule loss due to inelastic and chemical processes. 
Such dissipative processes are key to dissipative quan- 
tum phase transitions, which is a maior area of interest 
in quantum many-body theory (3Q|-[32|. 

We acknowledge useful discussions with Immanuel 
Bloch, John Bohn, Silke Ospelkaus, Luis Santos, and Pe- 
ter Zoller. This work was supported by the National 
Science Foundation under Grant PHY-0903457. 



Appendix A: The Internal Hamiltonian 

A ^ E molecule in its electronic and vibrational ground 
states has three angular momentum degrees of freedom: 
the rotational angular momentum N and the nuclear 
spins Ii and I2. In this work we shall use the coupling 
schemes | (hh) INFMf) and \I1M1I2M2NMN), which 
we refer to as the coupled and uncoupled bases, respec- 
tively. Explicit expressions for all single-molecule matrix 
elements in both bases are provided in Appendix O The 
relevant Hamiltonian for the internal degrees of freedom 



iJin may be written as a sum of rotational and hyperfine 
terms as 



where 

2 



(Al) 



(A2) 



Hm = ^CiN ■ I, + C3I1 ■ T ■ I2 + C4I1 ■ I2 + ^ V, ■ . 



i=l 



i=l 



(A3) 



The rotational term Eq. (|A2p corresponds to the Hamil- 
tonian of a rigid spherical rotor with (2A^ + l)-fold de- 
generate eigenstates \NMn)-, Mn being the projection of 
N on a space-fixed quantization axis [33|. The eigenen- 
ergies are given by Enmn = BnN {N where Bn is 
the rotational constant of the molecule (we use the nota- 
tion Bn instead of the more common B to avoid confu- 
sion with the magnetic field niagnitude B). In the case 
of '^^K^'^Rb, Bn=I.IU GHz The rotational level 

splitting defines the largest intrinsic energy scale for 
molecules. 

The first term of the hyperfine Hamiltonian, Yl'i=i ^^i^' 
Ii represents the interaction of the nuclear spins with the 
magnetic field created by the rotation of the molecule, 
and is governed by two coupling constants ck and CRb 
related to the nuclear shielding tensor. For ^^K^^Rb, 
these have been determined from density functional cal- 
culations to be ~20Hz and ~100Hz, respectively [23] • 
Because of the smallness of these constants and the fact 
that this term does not couple states with different 
this term plays a very small role in the spectra. 

The two nuclear spins have nuclear magnetic moments 
which interact via a resonant dipole-dipole interaction 



dd ^ 



Ii-l2 3(Ii-R)(R-l2) 



(A4) 



where gn is the proton (/-factor and R the vector joining 
the two nuclei [33]. This may be written as the contrac- 
tion of two rank-2 spherical tensors as 



dd = -g%^i% (Mo/47r) {R-')^/6{c f^ ■ (T {Iuh) f^ 

(A5) 



where (C)^ ^ is an unnormalized spherical harmonic in 
the relative degrees of freedom. The nuclear spins can 
also interact indirectly through the electron spins, and 
do so even for configurations [33|. This indirect in- 
teraction is represented by a tensor J which may be de- 
composed into its isotropic part Jiso and its anisotropy 
AJ = J\\ — J±. The combination of direct and indirect 
nuclear spin-nuclear spin interaction may thus be written 
as the sum of a scalar interaction and a tensor interaction 



^hf- 



dd ■ 



" -f^indi] 



C4I1 • I2 + C3I1 • T • I2 



(A6) 
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where C4 = Jiso, C3 = gjj/J^% {/iq/Att) {R ^) — A J/3, and where is a unit vector along the space fixed spherical 



the tensor T contains the angular dependence of the ten- 
sor interaction. C3 is of order lOHz for the various iso- 
topes of KRb, and so plays a very small role in the spec- 
tra. C4 splits the various levels according to their total 
nuclear spin / as 



{{hl2)INFMF\c^h ■ I2I {hh) I'N'F'M'p) 

= ^IJ'^N,N'^F,F'^MfM'f 



(A7) 



C4 is of order lOOHz-lOkHz for isotopes of KRb, and so 
is the dominant hyperfine contribution for TV = in 
the absence of external fields, see Fig. [71 Note that C4 
may be either positive or negative. For ^^K^^Rb, C4=- 
20.304kHz [23], and so the lowest energy states for TV = 
in zero field are the highest nuclear spin states / = 11/2. 

The final term in the hyperfine Hamiltonian is the in- 
teraction of the quadrupole moment of the nuclei with 
the gradient of the electric field produced by the elec- 
trons. We may represent this interaction by the sum 
Yld^i ' where is a second rank spherical tensor 
describing the electric field gradient at the i^^ nucleus 
and Qi is a second rank spherical tensor describing the 
nuclear quadrupole of the i^^ nucleus. The pertinent cou- 
pling constants (eqQ)- which arise in the matrix elements 
of this Hamiltonian are of order 100-lOOOkHz, making 
it the largest term in the hyperfine Hamiltonian. The 
quadrupole term doesn't affect the TV = level, however, 
and so the scalar spin-spin coupling dominates there. In a 
strong DC field the rotational levels become deeply mixed 
and the nuclear quadrupole thus becomes the dominant 
hyperfine contribution for all states. 



Appendix B: Interactions with static external fields 

Polar molecules such as heteronuclear dimers can cou- 
ple to external fields either through their permanent elec- 
tric dipole moment, through magnetic moments gener- 
ated from their rotation or nuclear spin, or through their 
polar izability tensor. The Hamiltonian representing in- 
teraction of the molecule with a static DC electric field 
Edc and a static magnetic field B may be written 

2 

i^F = - d • Edc - ^r^ArN • B - ^ gi/J^Nh • B (1 - di) . 

i=l 

(Bl) 

For molecules the permanent dipole moment d lies 
along the internuclear axis which defines the p = axis in 
a spherical coordinate system rotating with the molecule. 
Because this basis leads to anomalous commutation re- 



lations [Ji^Jk 



-ihcijkJk [3J we find it convenient to 



transform to the space-fixed frame where the angular mo- 
mentum operators satisfy the normal commutation rela- 
tions [Ji, Jk] = ihcijkJk, giving d • = = dC^^^ (6>, 0), 



q direction and Cc 



(1) 



is an unnormalized spherical 



harmonic whose arguments and (j) are the polar and az- 
imuthal angles of the internuclear axis in the space fixed 
frame. Taking matrix elements of dq in our two basis sets 
yields 

{hMihM2NMN\dq\hM[hM^N'M'j^) 
= SM,Mi^M,M^dy/{2N + 1) (2TV' + 1) (-1)^" 
f N 1 N'\ f TV 1 N' \ 







-Mn q M'j, 



(B2) 



N 



{{hh) INFMF\dq\ {hh) hN'F'M'p) 

^2F-Mf+I+N'+N+1 N 1 N' 



5ij^d{-l)' 







X V'(2TV + 1) (2TV' + 1) {2F + 1) {2F^ + 1) 

^ f TV lTV'\fTVF/l 
i -Mn q M'^ \ F' N' 1\ ' 



(B3) 



where 



Jl J2 J3 

mi 7712 rris 



is a Wigner 3-j coefficient and 



J Jl 32 J3 I is a Wigner 6- j coefficient [35^. We see that 
[ j4 js je J 

the rotational eigenstates have no net dipole moment, 
but that the dipole operator couples the state |TV, F, Mp) 
with the states |TV±1, F±l, Mp^q). The introduction of 
a DC electric field Edc with Hamiltonian — d • Edc cou- 
ples these levels and induces dipole moments, breaking 
the rotational symmetry and removing the (2TV + l)-fold 
degeneracy. Typical molecular dipole moments are mea- 
sured in Debye (D), where lD=503.4MHz/(kV/cm), and 
so the DC field becomes the dominant contribution to the 
Hamiltonian for modest fields of a few kV/cm. The per- 
manent dipole moment of KRb has been experimentally 
determined to be 0.566D [3]. 

On the scale of the rotational constant, the effect of a 
DC field on the single-molecule energy spectrum is as in 
Fig. [51 It is quadratic for field energies small compared to 
the rotational energy but becomes linear in stronger fields 
because the field strongly mixes states of opposite par- 
ity [sil . We consider the quantization axis to lie along the 
field direction, and so states with the same value of \Mn\ 
remain degenerate. A universal plot for all molecules 
results on this scale if the energy and field strength (i£^DC 
are both scaled to the rotational constant. 

The average orientation of the molecule with the elec- 
tric field can be obtained with the Feynman-Hellman the- 



dE 



d{dEBc) ' 



(B4) 



where E is the energy eigenvalue. The energy eigenvalue 
is dominated by the GHz scale structure, thus the degree 
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Electric Field (kV/cm) 

FIG. 5: GHz scale view of the Stark effect for KRh. Introduc- 
tion of a DC field breaks the degeneracy between all states 
with the same N but different | Mn \ • The large electric dipole 
moment causes GHz scale energy shifts which completely ob- 
scure the hyperfine splittings on the scale of this plot. Because 
the dipole moment is the same for any isotope of KRb, the 
Stark effect on this scale is the same for all isotopes. 




'•^0 10 20 30 40 50 

Electric Field fkV/cm) 



FIG. 6: Induced dipoles for KRb in an electric field. The 
N = and = 1, Mn = ±1 levels orient along the field, 
giving rise to positive dipole moments. The A/" = 1, Mn = 
state antialigns with the field for small fields, but aligns in 
stronger fields. All resonant dipole moments approach the 
"permanent" value 0.566D as the field strength increases. 



of alignment with the field is essentially independent of 
the hyperfine structure. From the degree of orientation 
we can also determine the effective space-fixed dipole mo- 
ment as d{cos 0). Fig.[6]shows the behavior of the induced 
dipoles as the field strength is increased. For all field 
strengths the N = {) and A^ = 1, Mn = ±1 states align 
with the field and so have a positive induced dipole mo- 
ment. In contrast, the A^ = 1, Mn = state antialigns 
with the field for weak fields and aligns with the field for 




10 20 30 40 50 

Electric Field (kV/cm) 



FIG. 7: kHz scale view of the Stark effect for ^^K^^Rb, 
N=0. All energies are shown relative to the GHz scale field- 
dependent average energy for A/" = 0, see Fig. [S] The inset 
shows the weak field region where the scalar spin-spin inter- 
action has split the levels according to / (equivalently, F), 
with larger / having lower energy. As the field is increased 
the nuclear quadrupole couplings split according to M/, and 
in large fields MRb and Mr also become well defined. See 
text for details. 



stronger fields. 

The magnitude of the field energy completely obscures 
the hyperfine splittings, and so to see the effects of hyper- 
fine structure we subtract from each state with a given A^ 
the field-dependent average energy of all hyperfine states 
with the same A^. For N = the results are shown in 
Fig. [71 For low fields the hyperfine splittings are dom- 
inated by the scalar spin-spin coupling and are of or- 
der C4, a few kHz. As the field is increased the vari- 
ous hyperfine states split according to |M/|. For large 
fields Ml and M2 also become well defined, which oc- 
curs because the energetic differences between states with 
AMn = ±1 become larger than the quadrupole coupling 
constants (see Eq. (|C9p ). Pairs of Mi and M2 which have 
the same |Mi + M2I are degenerate, and the state with 
I Ml + M2 1 = is degenerate due to reflection symmetry 
in the plane of the electric field vector. 

Because of the signs of the quadrupole couplings for 
^^K^'^Rb, the lowest energy states are those with MRb 
the largest and Mk the smallest. Because the kHz scale 
Stark effect depends on several molecular parameters it 
cannot be put into a universal form for all molecules 
like the GHz scale Stark effect. However, the qualitative 
structure will be similar for all molecules with nuclear 
quadrupole couplings; key differences being the number 
of nondegenerate levels and the energetic ordering of the 
magnetic quantum numbers [l^. The hyperfine Stark 
effect for A/" = 1 and other molecular species as well as 
the effects of electric fields on microwave spectra may be 
found in Ref. 0. 

Magnetic fields couple to the magnetic moments gen- 
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2000 




Magnetic Field (G) 




Magnetic Field (G) 



FIG. 8: Zeeman effect for K^^ Rb, N=0. The magnetic field 
splits the hyperfine levels according to their projections Mk 
and MRb with splittings between adjacent levels of order kHz 
for the experimentally relevant range B ^ 550G. The lowest 
(highest) energy state corresponds tom^ = —4 + 3/2 = —5/2 
(5/2). The zero field splitting is set by C4 and is not visible 
on the scale of this plot. 



FIG. 9: Zeeman effect for ^^K^^Rb, N=l. The zero field 
splitting is caused mainly by the nuclear quadrupole intrac- 
tion and separates the levels into groups of well defined F. 
The much larger zero field splitting causes avoided crossings 
between states with the same Mp to occur at much higher 
fields than in the N — case. 



erated by the rotation of the molecule and by the nuclear 
spins. The former interaction is given by —grfiN^ • B, 
where gr is the rotational g factor of the molecule and /i n 
is the nuclear magneton efi/2mp=762.259Hz/G [37]. The 
latter interaction is given by — ^i/^A^Ii ' B (1 — (7^), 
where gi and ^2 are the g-factors of nucleus 1 and 2, 
respectively, and is the isotropic part of the nuclear 
shielding tensor for nucleus i. The rotational contribu- 
tion is typically much smaller than the contributions from 
the nuclei, due to smaller g-factors and the fact that 
the isotropic parts of the nuclear shielding tensors are 
typically only a few parts per thousand. For example, 
in "^^K^^Rb gr = 0.0140, gK = -0.324^^/Rb = 1.834, 
ctk = 1321ppm, and <jRb = 3469ppm [23 • We neglect 
diamagnetic contributions to the Zeeman effect, as these 
contributions are small for the fields we consider. 

Typical experimental magnetic fields are ~ 550G be- 
cause of the Feshbach association stage of the STIRAP 
procedure [2^]. In Fig.[8]we show the Zeeman effect for the 
N = level of "^^K^^Rb for fields up to this range. We 
see that the magnetic field splits the spectrum accord- 
ing to the nuclear spin projections Mk and MRb, with 
larger (smaller) Mk (MRb) having lower energy due to 
the signs of the ^-factors for ^^K^^Rb. Because the nu- 
clear quadrupole interaction doesn't affect the = 
level, the zero field splittings are determined by the small 
scalar spin-spin coupling parameter C4. The Zeeman term 
dominates over the scalar spin-spin coupling at very low 
fields and so the effects of the scalar spin-spin coupling 
are not discernible on the scale of this plot. Additionally, 
the Zeeman contribution at these fields is larger than the 
hyperfine Stark splittings from the largest electric fields 



accessible in current experiments, see Fig. [71 

The spectrum for the A" = 1 level of ^^K^^Rb is shown 
in Fig. [HI It is greatly complicated by the fact that 
there are three times as many states as the A" = case 
(corresponding to the allowed Mat). Also, the nuclear 
quadrupole interaction affects the A" = 1 level, caus- 
ing the large zero field splittings. These larger zero field 
splittings delay the separation of the levels into well de- 
fined Ml and M2 , and also causes a complicated series of 
avoided crossings between states with the same Mf- 



Appendix C: Explicit values for the single-molecule 
matrix elements 

Here we present the matrix elements of the single- 
molecule terms of the Hamiltonian ()Aip in the coupled 
and uncoupled basis sets. We adopt the conventions of 
Zare [35|. 

The matrix elements of the rotational Hamiltonian are 
given by 

{IiMiI2M2NMn\BnN^\IiM[I2M^N'M'j^) 

BnN{N+1), (CI) 



{{hh) INFMF\BN^'^\{hl2) I' N' F' M'p) 

= Si,rSN,N'SF,F'SMF,M'^BNN {N + 1) . (C2) 

The matrix elements of the rotation-spin Hamiltonian 
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are given by 



are 



(/iMi/2M2A^Mjv| ^CiN • Ii\hM[l2M!,N'M'^) 



^q+N-Mr, I N 1 N 

I -Ma, q M'j^ 



1 h 

Mi -q M[ 



X ^/N {2N +1){N + 1) h {2Ii + 1) {li + 1) , (C3) 



{{hh) INFMf\ CiN • 1,1 (J1/2) I'N'F'M'p) 

i=l 

r / .-.I+N+F+h+h + l \ I N F \ 

= Sn,n' (-1) I ^ ^, ^1 

X x/A^ {2N +1){N+ 1) (/ + 1) {21 + 1) (2/' + 1) 



- ' ^ ^ V/i(2/i+l)(/i + l) 



^ /' /i 1 ^ 

+ <5m„Mi(-1/c2| I V/2 (2/2 + 1) (/2 + 1) 



(C4) 



The matrix elements of the scalar spin-spin coupling 
are 



{IiMiI2M2NMn\c31i ■ T ■ hlliM^hMl^N' M'j^) 
= -csVQ ( ^ q ^' ) Vm+1){2N' + 1) 



X ^/IlI2 (2Ji + 1) (2/2 + 1) ih + 1) ih + 1) 



q-M7V+/l-Mi+/2-M2 I N 2 N' 

\ -Mn q M'^ 



^(l,m;l,-g-m|2,-g) I 



-M2 -q-m I ' 



h 1 /i 



(C7) 



((/1J2) JA^FMFlcali • T • I2I (J1/2) I'N'F'M'p) 
.I'+F i I N f\ 



X V(2iV + 1) (2iV' + 1) 



[ A^' /' 2 

A/" 2 A/"' ^ 
0, 



/l /l 1 
h h 1 
I r 2 



X v/30 (2J+l)(27' + l)/i/2 

X v/(/i + 1) {h + 1) (2/1 + 1) (2/2 + 1) . (C8) 



{IiMj2M2NMN\ciIi ■ l2\IiM[l2M;,N'M'f^) 



I1-M1+I2-M2 



X VC^h + 1) h (/i + 1) (2/2 + 1) h {h + 1) 



' \-M^ q M[ 



h 1 h 

-M2 -q M'2 



(C5) 



((/1/2) /iVFMF|c4li • I2I (/1/2) I'N'F'M'p) 

= ^I,r^N,N'5F,F'^MF,MF 

X |^[7(7 + l)-7i(Ji + l)-J2(/2 + l)] . 



(C6) 



The matrix elements of the tensor spin-spin couphng 



The matrix elements of the nuclear quadrupole Hamil- 
tonian are given by 



(/iMi/2M2iVM^| ^ V, • Ql^\hM[hM'^N' M'^) 



2=1 q 

( N 2 N' \ ( h 2 Ii\ 



-Mn q M'j^ J ^ 
N 2 N'\ f L 2 1, 



Mi -q Ml 

1 







-li li 



(C9) 



13 



((/i/2)/7VFM^| ^ V, • Q,| {hh)I'N'F'M'p) 

i=l 

— of,f'OMf,Mf^ V -I-J 
X ^(2^" + 1) {2N' + 1) {2r + 1) (2/ + 1) 
' N 2 \ i I N f] 



X 



^ y [ 7VW 2 



/Wi 2 I V -h /i 

+ J,,,,,(e,Q),(-l/j 



-/2 lo 
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